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Abstract 

This paper discusses and compares six different methods for 
calculating the Jacobian for a general N-degree-of freedom 
manipulator. We enumerate the computational efficiency of 
each in terms of the total number of multiplications, addi- 
tions /subtract ions, and trigonometric functions required as 
well as in terms of the number of matrix-vector operations 
needed. We also give the execution times on a POP- 1 1/70 
minicomputer for determining the Jacobian for an example 
seven-degree-of freedom manipulator. This paper formulates 
one of the best new methods for determining the Jacobian. 

1. Introduction 

The control for most state-of-the-an industrial robots 
is relatively simple, based on servos at each joint. 
However, the controls for the next generation of robots 
must allow adaptation to the environment through the 
use of various kinds of sensory information (Dodd 
and Rosso! 1979), For this case, control is most natu- 
rally implernen ed in Cartesian coordinates at the 
end-effector. 

Some kind o transformation is needed between end 
effector coordirates and joint coordinates. An inverse 
kinematics procedure may be used to compute joint 
angles when the end-effector position and orientation 
are known. However, this may result in extremely 
complex equations that are difficult to derive and to 
implen;?nt in r?al time for control (Duffy 1980). Fur- 
thermore, for t le case of a redundant manipulator 
(more than 6 d grees of freedom), the inverse kinema- 
tics equations ? e in general underspecified and other 
approaches lo :ie conir-jl based on the Jacobian liave 
been proposed ,Liegeois. 1977: KJein and Huang 1983; 
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Trevelyan, Kovesi, and Ong. 1983; Yoshikawa 1983). 

The Jacobian relates joint rates to end-effector rates 
and may be used if control is based on resolved rate 
(Whitney 1969) and if position feedback is closed at 
the end-effector (Klein and Briggs 1980). Also, the 
transpose of the Jacobian relates end-effector forces to 
joint torques and may be used when force control is 
needed (Wu and Paul 1982). 

This paper compares several approaches that have 
been proposed for computing the Jacobian, including 
a new approach introduced here. These approaches 
are compared in terms of their computational effi- 
ciency since the Jacobian must be computed in real 
time for control An enumeration is made of the total 
number of multiplications, additions/subtractions, 
and trigonometric functions that are required to deter- 
mine the Jacobian. These are given as a function of 
the total number of degrees of freedom, N. The best 
approaches give computation that varies linearly with .V. 

The execution times for determining the Jacobian 
(using PASCAL on a PDP-l 1/70 minicomputer) for 
an example seven-degree-of-freedom manipulator is 
also presented. Finally, the number of matrix-vector 
operations involved in computing the Jacobian for 
each approach is enumerated. These may provide the 
most important basis of comparison for future matrix- 
vector processors. 



1.1. Notation 

The basic parameters used to describe the kinematics 
of a manipulator, as first presented by Denavit and 
Hailenberg ( 1955) and extended and given in detail in 
Paul (1981) are used throughout. In particular, a coor- 
dinate system is attached to each link of the manipula- 
tor with the z-axis directed along the joint axis. The 
links are numbered from 0 at the base to A" at the 
end-effector. A separate coordinate system, labeled 
with the subscript A* + I (or E for end -effector) is also 
fixed to the end-effector at any desired point. 
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Fii(. L Link parameters Fig. 2. Depiaion of vectors 

associated with link i. Pj> rj. pf, and transform 

'~^U, for a manipulator 




The relative position of joint /, i?,, is with respect to 
the z-axis of the previous link, z,_i . That is, links / — I 
and / are connected at joint /. 

Four parameters are used to. describe each successive 
joint and link pair— the joint angle (5,) and offset 
distance {d^) as well as the link length (a^) and twist 
(a,) (see Fig. 1 ). From these parameters, the 4X4 ho- 
mogeneous transformation, '~*T,-, which relates posi- 
tions in coordinate system / to those in coordinate 
system / — 1, may be computed. For cither a rotational 
or a sliding joint the result is as follows (Paul 198 1 ): 




X 



position vectors from the end-effector to link 7, r^, and 
from the base to link j, p^. These are all shown in Fig. 2. 



'T; = 



cdi —sdiCai sdiSUi 
sOj cdjCai —cOfSai 
0 5a, cai 



0 



0 



0 



1 



(1) 



where cO^ and 50, indicate the cosine and sine, respec- 
tively. The transformation between other link coordi- 
nate systems may be obtained through multiplication 
of the intermediate transformation matrices. 

The '~*T, matrix gives both position and orientation 
changes between successive coordinate systems. The 
lop-left 3X3 part of the matrix gives all orientation 
information and will be denoted as '"'U, . The first 
three elements of the right-hand column give all the 
relative position information and will be denoted as 
'~'pf . Also of value in describing the various ap- 
proaches to computing the Jacobian are the use of 



2. The Jacobian Matrix 

The firiarforrh of the Jacobian matrix, as determined 
in the various approaches, may differ in two impor- 
tant ways: ^ 

1. The components may be expresse'd in any 
coordinate system from 0 to iV 4- 1 . 

2. The reference point on the end-effector (either 
real or ficticious) for which the translational 
velocity is instantaneously computed may be 
chosen somewhat arbitrarily. 

These lead to the following notation for the Jacobian 
matrix, J, and the variables it relates: 



(2) 



where 
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% 



(3) 



and 0) is the rotational velocity vector of the end-effec- 
tor, V is the translational velocity vector of the end-ef- 
fecton and q is the vector of joint rates. The compo- 
nents of the end-effector velocity and the Jacobian are 
expressed in the ^th link coordinate system, as indi- 
cated by the leading superscript. The reference point 
for the end-effector velocity is understood to be in- 
stantaneously at the origin of coordinate system m, as 
indicated by the trailing subscript. For m equal to 
other than A* or ;V + I . the velocity reference point is a 
fictitious point not physically on the end-effector. 
However, it may be chosen in this way to reduce the 
computational complexity of determining the Jacobian. 

In the following sections, six different methods for 
computing.the Jacobian matrix are presented. In all 
cases, the fundamental approach from the appropriate 
reference is given. The notation, however, is standard- 
ized, and nonessential details (for purposes of this 
comparison) are eliminated. 

Several of the methods for computing the Jacobian 
are based on earlier work by Pieper (1968) and Whit- 
ney (1972). In fact, in many cases the latest work is 
just a detailed elaboration of the previous work. Hov\- 
ever, since the exact computational approach is not 
made explicit in Pieper and since closer attention to 
computational detail has resulted in somewhat greater 
efficiency than that presented in Whitney, these eariier 
methods will not be compared in this paper. 



For a rotational joint, the angular velocity of link / 
may be computed from that for link / — 1 and the rel- 
ative rate at joint /: 



The appropriate component equation is 



(5) 





"o" 






0 






_ 1_ 





The general form for ' *aj,-_, , 



may be used in Eq. (6), with the following result: 



(6) 



(7) 



(8) 



Equating the coefficients for Eqs. (4) and (8) results 
in a recursive method for computing % . Using a simi- 
lar approach for the translational velocity (with coeffi- 
cients 'Pj) and considering the case of sliding joints 
result in the following sets of recursive equations: 



0 
0 



(9) 



N) revolute joint: 



2. 1. Method I 

In the first method, as presented by Vukobraiovic a: J 
Potkonjak (19791 the angular velocity of link / is ^-.i it- 
ten as a linear function of the previous joint \ clociti(. s: 



/ = 1 . 2. 



- 1 1 : r I- 1 



7=1.2. 



= 1 



;V) sliding joint; 



. . / - 1, 
/V + I ; 



(10) 



(11) 



'(Oi=2'y/lj- 
y-i 



(4) 



'A =■= 'y, X 'pf . 



/=1,2 /V) revoluie joint; (12) 
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.7? = .-n_7 



0 
0 

1. 

/■= 1,2, 



(13) 



. yV) sliding joint; 



'y?, = '-'LT-')5f, + ('y;X'pr), 

;= 1, 2, ...,/- 1, 
/=2 A^-f- 1. 



(14) 



For / = A -f 1 in Eqs. (11) and ( 1 4), the coefficients 
as determined are just the elements of the Jacobian 
matrix, -^"^^Ja + i = ^Jr- This rnethod is rather ineffi- 
cient because for / == A+ 1, A+ 1 Jacobians are com- 
puted. First, the Jacobian is determined for a manipu- 
lator consisting of the first link only, then for a 
manipulator consisting of the first two links, then 
three links, and so forth up to links plus the end-ef- 
fecior. The recursions in all cases are from the base to 
the end-effector. This is a rather inefficient method 
because all that is needed is the Jacobian for the entire 
manipulator. It may be very useful, however, in the 
study of dynamics, where it is necessar\^ to compute 
the Jacobian for each set of links in the manipulator 
(Vukobratovic and Potkonjak 1979). 



2.2. Methods II, IlL and IV 

The next three methods use the same basic concepts 
but result in different forms for the Jacobian matrix. 
They are based on earlier methods by Pieper (1968) 
and Whitney (1972). With close attention to computa- 
tional detail, they provide relatively efficient methods 
for computing the Jacobian. 

The angular velocity for the end-effector may be 
written as a function of the joint rales (for rotational 
joints): 



(15) 



The translational velocity for the end-effector may be 
determined by noting that a particular component is 
just the appropriate cross-product of the individual 
joint rate vector and position vector from the joint axis 
to the velocity reference po\n\ (-t,). That is, 



v = -(qiZ,;)Xto-(q27.i)Xt,-- • • 

-(q.vZA-i)XtA-.* ^ • 

Comparing Eqs. (2) and (3) with Eqs. ( 1 5) and ( 1 6), 
we find that the coefficients of the Jacobian matrix, 
for the case of a rotational joint, may be written in 
terms of the joint axis vector and relative position vec- 
tor: 



A = -(z,_, Xt,_,) 



rotational joint. (17) 



For the case of a sliding joint, similar expressions may 
be derived: 



}', = 0. 
A = z,. 



sliding joint. 



(18) 



The three methods presented in this section differ in 
the coordinate system used to express the components 
of the Jacobian and also in the velocity reference point 
chosen. In Waldron's work (1982). the components of 
the Jacobian are determined with respect to the base 
coordinate system 0. and the velocity reference point is 
chosen to be the origin of the base coordinate system. 
°Jo is thus computed. Since p, is defined as the posi- 
tion vector from the origin of the base coordinate 
system to that of link / (see Fig. 2).. then the following 
recursive equations may be used: 



°Uo = I: 

/■= 1, 2. . 



,A^- 1; 



°Po = 0; 



0 
0 
1 

/"= 1,2, . . 
/= 1,2, . . 



, A') revolute joint; 



, A') sliding joint; 



(19) 



(20) 



(21) 



(22) 



(23) 
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v.=°p,_.+°u,'pr, 

/= 1,2, 

°A=°y.x(-°p.-i). 

/= 1,2, 



, A^- 1; 



(24) 



, N) revolute joint; 



0 
0 

1 



(26) 



, N) sliding joint. 



Note in the above equations that the components of 
z,_, in coordinate system / - I are [0 0 1]^ Also, I 
is the identity matrix. In Eqs. (20) and (24), the recur- 
sions for computing the orientation matrix and posi- 
tion vector are directed from the base out to the end- 
effector. 

In Olson and Ribble^s work (Ribble 1982), the com- 
ponents of the Jacobian are determined with respect 
to the base coordinate system, and the velocity refer- 
ence point is chosen to be the center of the end-effector 

^J^ is thus computed. Since r, is defined as the posi- 
tion vector from the origin of the end-effector coordi- 
nate system to that of link /, Eq. (25) may be replaced by 



0- = On _ On 

r,-i — Pi- 1 Pa^+i 



i= \, . . . ,N. 



(27) 
(28) 



Also, the indices on Eqs. (20) and (24) must be 
changed to A' + 1 . In Olson and Ribble's work the 
recursion for computing the orientation matrix is from 
the base out. The recursion for the position vector is 
also from the base out; however, Eq. (27) must be 
added to give the velocity reference point at the center 
of the end-effector. While this method is slower than 
that of Waldron (1982), an advantage is that the veloc- 
ity reference point is a real point on the end-effector 
and so the values of °A relate directly to the velocity of 
the cnd-cffector. 

In Renaud's work (1981), the components of the 
Jacobian and the velocity reference point are both 
associated with a link which is appro.ximately mid- 
way between the base and end-effector. 'Jf is thus 
computed. The recursive equations for computing y 
and are 



(29) 



z = ^+ 1,^ + 2, 



,N- 1; 



(25) 'U,_, 



i==e, (£-1), 



. 1; 



'y. = o. 



0 
0 

I 

i = r, 2, 

/ = 1,2, 



, N) revolute joint; 
, A''} sliding joint; 



0; 



/• = ^+ 1,^ + 2, . . . ,A^- 1; 



/ = ^,(^-l), 

/•= 1,2, . . . 



. . . , 1; 

, N) revolute joint; 



0 
0 
1 

/= f, 2, 



(30) 
(31) 

(32) 

(33) 
(34) 
(35) 

(36) 

(37) 

(38) 



, N] i liding joint. 



The position vector t, in the previous equations is 
from the origin of coordinate system m = ^ to the ori- 
gin of coordinate system /. Also note that ttie recur- 
sions are from link ( in toward the base and out 
toward the end-effector. 

The main advantage of this method is lh;U the ma- 
nipulator is broken into two manij ulators of approxi- 
mate length N/2, and the Jacobian 'or a manipulalor 
of length N/2 is calculated twice..! le matripes at the 
beginning of the Jacobian calculat' .)n are silnple. and 
Renaud's method, starting twice, t kes adve ntage of 



\ 



2.3. Method V 



The method described in this section is ricv. , rcsuitini^ 
from an eftbrt to compute in the most efficient 
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\va\- possible. We see from Eqs. (9) through (14) that 
■ Method I is somewhat inefficient because the elements 
of ;V different Jacobiansare computed. That is, the 
Jacobian for a manipulator consisting of the first hnk 
only is determined, then for a manipulator consisting 
of the first two links, then the first three links, and so 
on. Also, the recursions are from the base to the end- 
effector. 

If the recursion is redirected to be from the end-ef- 
fector to the base and if values of y and not needed 
are not computed, then the resulting equations are 



A'+ 1 1 1 = 1 • 



,2, 1: 



;V-rlv =: 



/ = 1 . 2 



A'') revolute joint; 



y, = o, 



/= 1,2. 



, A') sliding joint; 



■ ,V+ 1 



= 0: 



A- I, 



/ = (A-f 1), ... ,2, 1; 

,v+iy_.X(-A+'r,_,), 
; = 1,2, . . . , A'} revolute joint; 



(39) 
(40) 

(41) 

(42) 
(43.) 
(44 J 

(45) 
(46) 



/= 1, 2, . . . , A'} sliding joint. 



It should be noted that the resulting equations are 
similar to Waldron's (1982) except that the recursion is 
in the opposite direction. This was not immediately 
apparent during the development of the equations. 
Once the common notation had been applied to all the 
methods, we recognized the similarity. The velocity 
reference point is now a real point on the end-effector, 
and the components of the Jacobian are expressed in 
the end-effector coordinate system. 



2.4. Method VI 

The final method comes from PauFj work ( 1 98 1 ) and 
also computes '-"J^, Paul, however, u ;es the 4 X 4 ho- 
mogeneous transformation matrix aj the basis for the 
computations The appropriate equa:ions arc 



J- IT = /—IX 'T 

* A + I * / * .V-»- 1 > 

'■ = (A'+ 1). 



2. 1, 



i-iijr 



/= 1. 2. 



.V+l,_ 



y, = 0, 



/= 1. 2. 



,N) revolute joint: 



A ) sliding joint: 



;=i.2,3. 

/ = 1. 2. . . 



0 
0 

1. 

revolute joint: 



(47) 
(48) 

(49) 

(50) 
(51) 



i-lT!?- 

/= 1. 2 



(52) 



, A'} sliding joint. 



Note that '"*U:{v.+ , is just column j of the 3 X 3 orien- 
tation part of '"*Tj^.+ , . Also note that^^^^r/,, is just the 
first three elements of the fourth colcfrtin of '"*T;^r+( . 
Furthermore, is just the j\h component of the 
vector ^"^'^y. 

It may be noted that the resuhs for Methods V and 
VI are quite similar except that in Method VI, the 
3X3 orientation matrix and 3 X 1 position vector of 
Method V are grouped into the 4 X 4 homogeneous 
transformation matrix form. If the bottom row of the 
homogeneous transformation matrix is understood to 
be composed of three O's and a 1, this seems to make 
little if any difference in the overall computational 
complexity. However, another difference is that the 
ordering of matrix multiplication in Eq. (48) is re- 
versed from that of Eq. (40), and this does result in 
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Table 1. Number of Computations for Each Method Maintaining the Natural Jacobian Form 
(A^ = Number of degrees of freedom) 









Addition/ 


Sines/ 




6 


Method 


J 


Multiplication 


Subtraction 


Cosines 


M 




Vukobratovic/Potkonjak 




iON^- \5N-\-9 




IN 


< 


lb 


/•n Waldron 


% 


30A^-55 


i5yv-38 


IN -2 






. HI Olson/Ribble 




30A^- 11 


18A^-20 


IN 






IV Renaud 


'J/ 


30N~S1 


\5N-66 


2N-2 






V Orin/Schrader 




30/^-25 


\SN-1S 


IN 






VI Paul 




30N- 18 


\AN- 15 


IN 


'% '2- 





i3 
/o 



2M 



increased complexity in Eq, (5 1) as opposed to that of 
Eq, (45). It appears then that some computational 
savings may be made in Paul's method if the order of 
matrix multiplication in Eq. (48) is changed. 

3. Comparison of the Methods 

We now assess the numerical efficiency of computing 
the Jacobian according to each of the six methods just 
described. In particular^ we enumerate the lolal num- 
ber of additions/subtractions, multiplications, and 
trigonometric functions performed in the process of 
evaluating the Jacobian. The results, given as a func- 
tion of the total number of degrees of freedom A'', are 
shown in Table 1. 

The results in Table 1 were generated through the 
use of a program that considered the components of 
matrices and vectors to be composed entirely of four 
symbols. The symbol x indicated a nonzero value 
other than ± 1; + indicated a value of + 1, — indicated 
a value of - 1, and 0 indicated a value of zero. The 
program produced a symbolic Jacobian while counting 
the number of multiplications, additions/subtractions, 
and trigonometric functions used to produce the Jaco- 
bian. Additions, subtractions, and multiplications that 
included zero were not counted. Multiplications that 
included ± 1 were totaled separately but have been 
added into the multiplication formulas of Table 1. It 
has also been assumed that the twist angle of each link 
is 0' or ±90" since most manipulators have this type 
uf kinematic configuration. 

From Table 1, we see that the best method is that of 
Renaud (1981). It should be understood, however. 



Table 2. Execution Time for Each Method on the 
PDP-11/70 



Method 


J 


Execution Time 
(ms) 


1 Vukobratovic/Potkonjak 




23.0 


II Waldron 


'Jo 


12.8 


III Olson/Ribble 




16.2 


IV Renaud 


'J. 


12.9 


V Orin/Schrader 




15.5 


VI Paul 




16.9 



that to use this form of the Jacobian for resolved rate 
control it is necessary to associate the angular and 
translational velocity of the end-effector with coordi- 
nate system where £ - (;V 4- l)/2. 

Methods I, V, and VI all result in the same form for 
the Jacobian, ^J^. Of these, the new method proposed 
in this paper requires the least amount of computation. 

All of the methods discussed in this paper have been 
used to simulate a seven-degree-of-freedom manipula- 
tor that has a somewhat general kinematic configura- 
tion (nonzero joint offsets and link lengths). A 
PDP-l 1/70 minicomputer with floating-point hard- 
ware and the PASCAL language were used, and the re- 
sults are given in Table 2. The results show a close 
correlation with those of Table I. 

In Table 2, any changes in ordering from that antic- 
ipated by reviewing Table 1 probably result From dif- 
ferences in the amount of overhead involved in in- 
dexed addressing. Also, in evaluating the execution 
times given in Table 2, it should be understood that 
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Table 3. Number of Computations for Each Metihod IncI iding Conver^^ion to 
(/V=^ Number of degrees of freedom) 







Addition/ 


Sines/ 


Method 


Mullipacation 


Subtraction 


Cosines 


I Vukobratovic/Poikonjak 


IOA'2- ISA' + 9 


A-*2 + 7A' - 2 


2A* 


II Waldron 


54A'- 31 


33A'- 37 


2A- • 


III Olson/Ribble 


48A'-23 


30A^ - 32 


2A^ 


, IV Renaud 


54A' - 8 1 


30A' - 79 


2A' 


V Orin/Schrader 


30A - 25 


ISA— 25 


2 A' 


Vi Paul 


30A - 18 


14A*- 15 


2 A' 






13^1 ' 


7^ 



Table 4. Execution Time for Each Method on the 
PDP- 11/70 

(All Forms Further Converted to ^Jf). 



Method 



Execution Time 
(ms) 



I Vukobratovic/Potkonjak 

II Waldron 

III Olson/Ribble 

IV Renaud 

V Orin/Schrader 

VI Paul 



23.0 
20.1 
19.1 
20.1 
15.5 
16.9 



the absolute magnitudes have little significance. That 
is. little effort was made to reduce.these times. Instead, 
much' effort was made to*ensure '^fairness" so that - 
relative magnitudes are significant. 

Perhaps a more appropriate way to compare the 
various methods is to transform the Jacobian produced 
by each method to a common form. In this case, ^J^ 
will be used as the final form. To move the velocity 
reference point from the origin of coordinate system 
m to the end-effector, the following transformation is 
appropriate: 



I 



' m 



0 



' m 



(53) 



To conven the components of 'y.^ and '^,£ from the 
^th coordinate system to the end-effector coordinate 
system, the following transformation may be used: 



f\ rr 



0 



0 



£ J 



(54) 



The total transformation is just a projective transfor- 
mation of screw coordinates (Sugimoto and Matsu- 
moto 1983) and is given as follows: 



0 



0 




{,■>■ 








0 


' m 


! ^ 






' m 


0 







(55) 



The results for each method after conversion to ^J^ 
are given in Tables 3 and 4. From these tables, we see 
that the last two methods are the most efficient, with 
the one proposed in this paper being the best. 

It should also be understood that the results pre- 
sented in these tables are not completely definitive. 
Each method has its own natural output form for the 
Jacobian. and to depart from this to transform it to a 
common form biases the results toward those that 
naturally produce this common form. However, 
whether the natural or common form of the Jacobian 
is used. Methods V and VI are especially computa- 
tionally efficient. 
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Table 5. Number of Matrix- Vector Computations for Each Method 









Matrix- 


Matrix 


Vector 


Matrix 








Vector 


Product 


Addition/ 


Product 


Method 


J 


Cross-Product 


E*roduct 


(3X3) 


Subtraction 


(4X4) 


I Vukobraiovic/Pptkonjak 




l/2(A^ ^ 3;vr) 


N^ + N 




1/2(A'2 + A^) 




11 Waldron / 






N- 1 


A^-2 


A^-2 




in Olson/Ribble 






A^+ I 


A^ 


2N 




IV Renaud 






A^-2 


A^-3 


;V-3 




V Orin/Schrader 






A^ 


A^ 


/V 




VI Paul 




IN 











One final comparison among the methods is given 
in Table 5. In this case, the basic arithmetic operations 
are for vectors and matrices. This is appropriate be- 
cause processors of the future should be equipped with 
such instructions. The results are similar to those ob- 
tained previously, in that Renaud's method is the best. 
Of those thai result in ^Jf, the method proposed in 
this paper is the best, (A 4 X 4 matrix product of ho- 
mogeneous transforms is basically equivalent to a 
3X3 matrix product along with a matrix-vector prod- 
uct and vector add.) 

4. Summary 

Six different methods for computing the Jacobian 
have been presented and compared for computational 
efficiency. The results indicate that the new method 
presented in this paper is the most computationally ef- 
ficient when the Jacobian is based on end-effector 
coordinates. A method proposed by Renaud (1981) 
gives the best results when an arbitrary coordinate 
system is used. 

To use a Jacobian of any particular form for re- 
solved rate control or force control the end-eflfector 
rates or forces must be associated with the same coor- 
dinate system. The total amount ofcomputation re- 
quired often depends on the specifics of the application 
for which the Jacobian is used. At a minimum, though, 
the results of this, paper should give a general indica- 
tion of the total computation required. 

The results presented here assume the use of a serial 
processor. A pipelined linear array of processors may 
also be used but requires a new formulation of the 



problem. This is the subject of present investigations 
at The Ohio State University (Schrader 1983), and the 
present paper provides the foundation for this work. 
Advances in both algorithms and architectures prom- 
ise to allow a control based on the Jacobian to be 
more viable for real-time implementation in the future. 
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